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EFFECTS OF WING-FUSELAGE FLOW FIELDS ON MISSILE 
LOADS AT SUBSONIC SPEEDS 
By William J. Alford, Jr. 

SIMMARY 

The flow-field characteristics around a swept-wing airplane model at 
low subsonic speed are described, and the loads induced on a typical mis- 
sile model while operating within these flow fields are presented. In 
addition, theoretical flow fields are compared with experiment and are 
used in first-order estimations of the resulting Induced missile loads. 


INTRODUCTION 


The loading problems associated with the external storage of teinks, 
bombs, nacelles, £ind missiles have become increasingly serious with 
increase in airspeed, and knowledge of these loads is desirable in the 
design of the store-supporting members. A study of available data from 
wind-tunnel tests (refs. 1 to 5)> in conjxmction with theory (refs. 4 
and 5)^ has indicated that these loads are primarily due to the nonuni- 
form flow field generated by the airplane. 

The purpose of this paper is to describe these flow fields at low 
subsonic speeds and to indicate the magnitude of missile loads that may 
result becavise of the nature of these flow fields. In addition, compeiri- 
sons of theoretical flow fields with experiment are presented. These 
flow fields, both theoretical and experimental, are \ised in first-order 
estimations of the resulting missile loads. 


SYMBOLS 


a 


a 


I ■ 


e 


angle of attack, deg 

local angle of attack between xy-plane and local flow direc- 
tion, a - e, deg (see ref. 1) 

angle of downwash, deg 

angle of sideslip 
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local angle of sideslip between xz-plane and local flow direc- 
tion, 3 + a, deg (see ref. 1) 

angle of sidewash, deg (see ref. 1) 

angle of sweepback, deg 

Mach nvmber 

aspect ratio 

taper ratio 

local wing chord, ft 

longitudinal distance, ft 

lateral distance, ft 

vertical distance, ft 

free -stream velocity, ft/sec 

free-stream dynamic pressure, Ib/sq ft 
local dynamic pressure, Ib/sq. ft 

missile length, ft 
wing span, ft 

longitudinal location of chordwise vortices (fig. 8), ft 
circulation, sq ft /sec 
chordwise vortex index 
spanwise vortex index 

longitudinal perturbation velocity parallel to free-stream 
direction (fig. 8), ft/sec 

perturbation velocity normal to spanwise line of constant 
sweep, ft/sec 

lateral perturbation velocity normal to free-stream direction 
(fig. 8), ft/sec 
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Subscripts : 

c/4 at quarter chord 

c/2 at midchord 

CG missile center -of -gravity position 

TEST CONFIGURATIONS 


Presented in figure 1 are config\irations used in investigations, made 
at low speed in the Langley 500 MPH 7 - 10-foot tunnel, in which flow- 
field surveys were made \ising the survey rake shown in the center of the 
figvire. This rake consisted of six probes, each of which was instrxanented 
to measure local pitch and sideslip angularities and dynamic pressure. 
These surveys were made under the fuselage and at each One-quarter semi— 
g pan location of the swept -wing— fuselaige combination and at the one-half 
semispan location of the vinswept -wing-fuselage combination. The range of 
vertical, spanwise, and chordwise position investigated is indicated by 
the dots, which represent the leading edge of the survey rake. 

Static force and moment measurements were obtained on a typical mis- 
sile model, shown in the center of figure 2, for various locations within 
the wing-fuselage flow fields. The range of missile center-of -gravity 
location investigated is shown by symbols at the data points. 

The results of these investigations indicated that the induced longi- 
tudinal characteristics were, in general, qualitatively similar for both 
the swept and unswept conf igvirations , whereas the induced lateral charac- 
teristics were considerably more seriously affected for the swept-wing 
configuration. In view of this result, the discussion deals hereinafter 
with the swept configuration. 


MEASURED FLOW FIELDS 


In defining a flow field, or velocity distribution, it is necessary 
to know both the direction and magnitude of the local velocity vectors . 

For convenience, the flow directions are expressed as vertical and lat- 
eral inclinations and the magnitudes are expressed as dynamic pressure. 

The local angles of attack as measured about the one -half semispan 
location of the swept-wing— ftiselage combination are presented in figure 5* 
The conditions depicted, which are for an angle of attack of 8° and zero 
sideslip, axe equivalent to the 6g service load limit of a typical fighter 
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airplane for a velocity of 450 mph at an altitude of 20,000 feet. These 
data are presented in contour formj that is, as lines of constant an^- 
larity for the physical space siirrounding the airfoil profile. The mis- 
sile outline is superimposed on the contours to show its size relative to 
the local wing chord and to indicate the angiilarity gradients to which 
it is subjected. Regions where the local angles are greater than 8° indi- 
cate an upflow and regions where these angles are less than 8° are indica- 
tive of downflow. The angularity variation along the center line of the 
missile, which in this case was 15 percent of the local wing chord below 
the wing chord plane, is presented in the lower portion of figure 5. 

Two missile center-of -gravity locations are shown to indicate both the 
angularity magnitudes and the angularity gradient that exist along the 
missile center line. For the chordwise range indicated, the local angles 
of attack vary from 5° to 12°, a gradient of 7° along the missile center 
line extended. 

Presented in figure 4, for the same vertical plane, are the local 
sideslip angles. It should be noted that the perturbation velocities 
which generate these lateral angularities are in a plane normal to the 
plane of the figure. The local sideslip angles below the wing chord 
plane represent an outflow direction (toward the wing tip), and the 
angles above the chord ^lane indicate an inflow (toward the plane of 
symmetry). The sideslip angles that exist along the missile center line 
are presented in the lower graph and indicate a gradient of 6° along the 
center line extended. Comparison of these sideslip gradients with angle - 
of -attack gradients of figure 3 shows that they are of the same magnitude 
and could be considerably more important since the supporting pylon would 
have its least structural strength in the lateral direction. 

The local dynamic pressures referenced to free-stream conditions are 
presented in figure 5* The effects of the induced longitudinal velocities 
are as would be expected, since the dynamic pressures are increased above 
and decreased below the wing chord plane. Sizable gradients are again in 
evidence over the length of the missile center line extended. These data 
also indicate that, for positive angles of attack, the dynamic pressure 
can effect load reductions but, for negative angles of attack, would cause 
large load increases. 


MEASURED MISSILE LOADS 


The preceding discussion has attempted to define and illustrate the 
flow phenomenon existing around the airplane. Having shown that the mis- 
sile is forced to operate in regions of adverse flow, it would be desir- 
able to correlate the flow characteristics with the loads induced on the 
missile. 
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The total missile normal forces, in pounds, for a free-stream dynamic 
pressijre of kk^ pounds per sq^xare foot are presented in figure 6. These 
data were obtained by assvnning a full-scale missile to be stored externally 
on a typical fighter airplane at a velocity of 450 mph at an altitude of 
20,000 feet. The effects of compressibility are absent since the missile 
loads data were obtained at low speed. The normal-force load center , in 
percent of the missile length from the missile nose, is also presented. 
These parameters are given for various locations of the missile center of 
gravity in fractions of the local wing chord, relative to the leading edge 
of the local wing chord. For comparison, the isolated-missile character- 
istics, at an identical attitude, are shown by the dashed lines. The two 
missile center-of -gravity locations (fig. 6) indicated by the solid sym- 
bols are identical to those shown in the lower graph of figiare 3. It 
should be noted that the total normal forces are a resvilt of an integra- 
tion of the angularities (fig. 5) an<i dynamic pressiures (fig. 5) over the 
length of the missile. Considering the rearward location of the missile 
center of gravity (x/c » 0.45), "the missile is seen to be operating in a 
region of reduced angularity and consequently experiences a reduced normal 
force when compared to the isolated missile. Further examination shows 
that the tail is operating in a slightly higher angular region than the 
missile wing (fig. 3) and consequently the load center is drawn aft. When 
the missile center of gravity is in the forward position (x/c « -0.75), 
the missile is experiencing higher angiilarities (fig. 5) and. the resulting 
normal forces are considerably increased (fig. 6). For the range of mis- 
sile center-of -gravity location Investigated, the normal forces are ini- 
tially decreased approximately 30 percent and later Increased about 50 per- 
cent with a total load-center travel of approximately 10 percent of the 
missile length. 

Presented in figure 7 are the total missile side forces and side- 
force load centers as a function of missile center -of -gravity location. 
Recalling that the airplane -missile combination is at zero sideslip, the 
comparable isolated-missile characteristics would be zero. In order to 
show the relative magnitude of these side loads, the isolated-missile 
characteristics at 6 ° sideslip are shown as the dashed lines. Once again 
the solid symbols represent the missile center -of -gravity locations shown 
in the lower graph of figure 4. Since the missile was below the wing, the 
local sideslip angles are in an outward direction (fig. 4) and cause neg- 
ative side forces (force directed toward the wing tip, fig. 7)* The side 
loads and load centers are as would be expected since, when the missile 
wings are in the highest angular region (fig. 4, x/c sa O. 5 O), the load 
center moves forward (fig. 7 ) and when the tail is in the highest angular 
region (x/c fsi '_o. 75 ) the load center moves rearward with a total load- 
center travel of 50 percent missile length for the chordwise rsuige indi- 
cated. The region of greatest side load is encountered when the missile 
wings are in the highest angular region. 
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THEORETICAL FLOW FIELDS 


Once a problem has been defined and shown to exist, it then becomes 
desirable to have, or to formulate, procedures by which the individual 
components of the system can be studied. 

In the present instance, the chief difficulty appears to be in the 
severe singular gradients that are generated by the airplane. The niathe- 
matical models used to calculate these flow fields assumed a simulated 
lifting siurface (fig. 8), the appropriate airfoil -section singularity 
distribution (by method of reference 6), and simple sweep theory (fig. 8). 
The simulated lifting surface, shown in figure 8, approximated both the 
spanwise and chordwise distributions of vorticity by discrete horseshoe 
vortices. The spanwise vorticity distribution was represented by 10 horse- 
shoe vortices and the chordwise distribution was represented by four vor- 
tices of equal strength, the chordwise locations of which were determined 
from the familiar two-dimensional circulation distribution shown in fig- 
ure 8. The nonlifting (thickness) effects were determined from the 
source-sink distribution that satisfied the two-dimensional airfoil- 
section boundary conditions (ref. 6) and were corrected by simple sweep 
theory (fig. 8) to account for wing sweep. 

Typical results of these calculations compared with experiment are 
presented in figiares 9 and 10. Figure 9 presents the local angles of 
attack and figinre 10 presents the local angles of sideslip for the one- 

half semispan location =0.5^ of the swept -wing— fuselage combination. 

The experimental data are shown as the symbols and the theory is shown 
as the solid-line curves. As is seen, the agreement is poor for a = -8°. 
This disagreement is presumed to be due to the fact that the flow on the 
suction side of the airfoil assumes characteristics that are nonpotential. 
The agreement is good between theory and experiment (figs. 9 and 10) for 
a = 0° eind for all positive angles of attack up to a = 24°, where theory 
then overestimates the local angles of attack. This is rather surprising 
from consideration of the nonpotential nature of the flow on the suction 
side of the wing surface. Calculations made for the three-quarter seml- 
span location have shown that theoretical values obtained by using the 
theoretical span-load distribution overestimate the experimental values. 


ESTIMATED MISSILE LOADS 


With the flow-field characteristics known, the next step is to use 
them, in conjunction with the missile component characteristics, to esti- 
mate the airplane induced missile loads. This estimate has been made and 
figures 11 and 12 present sample comparisons of theory and experiment. 
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Presented in figure 11 are the estimated normal forces and normal- 
force load centers. Experimental loads are again represented by the sym- 
bols. The estimation obtained by using the experimental flow fields is 
shown as the solid curves and the estimation by using the theoretical 
flow fields is shown as the dashed curves. Good agreement in estimating 
the normal forces is obtained by both estimates over the rear portion of 
the chord, with evidence that theory gives values too low ahead of the 
leading edge. The load centers are also seen to be well predicted. 

In the case of the estimated side forces (fig. 12) both estimates 
are low, although the side -force load centers are well predicted. The 
reason for the lack of agreement for the side forces is not completely 
understood. 

The data that have been presented were obtained at low speed in 
order to permit convenient examination of the nature of the complex flow 
that exists around airplanes. The xise of a missile as the store configu- 
ration was for illustrative purposes, and the approach utilized should be 
equally valid for other external stores, although such application has 
not, at present, been adequately demonstrated. 

The results of a brief theoretical study have indicated that the 
effects of compressibility, for subcritical Mach numbers, are to generate 
larger flow distortions and consequently larger induced missile loads, 
although the flow structure remains similar to that of an incompressible 
nature . 


CONCLUDING REMARKS 


In summary, these results have shown that the flow fields in regions 
beneath the one -half semispan location of a swept-wing airplane model can 
be calculated. These flow fields can then be used in first-order estima- 
tions of the loads experienced by a missile while operating in these flow 
fields . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., April 25, I 955 . 
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FLOW-FIELD-SURVEY CONFIGURATIONS 

a=-8“T0 24»; M=O.I3 



Figxare 1 


MISSILE-LOADS CONFIGURATIONS 
a = -8“ TO 20*; M=0.I3 
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LONGITUDINAL- ANGULARITY FIELD 
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Figure 3 


LATERAL- ANGULARITY FIELD 
a = 8®; A=0® 




Figure k 
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SIDE FORCES 


a»8*i j9«0*i qo = 445 LB/SQFT 




Figure 7 


THEORETICAL FLOW-FIELD MODELS 

LIFTING SURFACE THICKNESS-SWEEP 




Figure 8 
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PREDICTION OF LONGITUDINAL ANGULARITIES 
■^■0.5; Z/C.-O.I5 



Figure 9 


PREDICTION OF LATERAL ANGULARITIES 

= 0.5i 2/c ■= - 0.i5j fi ‘0" 



Figure 10 
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ESTIMATED NORMAL FORCES 

qQ = 445 LB/SQFT 
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Figure 11 


ESTIMATED SIDE FORCES 

a-8*i A“0®j qQ«445 LB/SQFT 
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